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Near-field heat engines are devices that convert the evanescent thermal field supported by a 
primary source into usable mechanical energy. By analyzing the thermodynamic performance of 
three-body near-field heat engines, we demonstrate that the power they supply can be substantially 
larger than that of two-body systems, showing their strong potential for energy harvesting. Theo¬ 
retical limits for energy and entropy fluxes in three-body systems are discussed and compared with 
their corresponding two-body counterparts. Such considerations confirm that the thermodynamic 
availability in energy-conversion processes driven by three-body photon tunneling can exceed the 
thermodynamic availability in two-body systems. 

PACS numbers: 44.40.+a, 84.60.-h, 05.70.-a, 78.67.-n 


A heat engine may, in general, be conceived as a device 
that converts part of the heat coming from a hot source of 
energy into mechanical work throughout an appropriate 
conversion system [I]. In contactless devices this heat is 
transferred to the converter by radiation, only. At long 
separation distances the maximum power which can be 
transmitted is bounded by the blackbody limit [2]. On 
the contrary, at separation distances smaller than the 
thermal wavelength, heat can be transferred to the con¬ 
verter also by photon tunneling, so that the flux can be¬ 
come several orders of magnitude larger than in the far- 
held regim e, as shown theoretically [3J and experimen¬ 
tally j5MlO|. Furthermore, it could be shown experimen¬ 
tally and theoretically that near-held thermophotovoltaic 
conversion devices can be used to harvest this energy by 
transferring it towards a p-n junction 11, (T^]. Thanks 


to the tunneling of surface phonon polaritons (SPPs) 
supported by the primary source, this energy transfer 
is quasimonochromatic, which is very advantageous for 
the energy conversion with a photovoltaic cell 12]. How¬ 
ever, recent theoretical works [Tdl | have demonstrated the 
strong potential of near-held heat exchanges in the most 
general context of heat engines for capturing low-grade 
waste heat for power generation. Beside this result, a 
mechanism of photon-tunneling enhancement has been 
predicted in three-body (3B) systems 14-1?} when pas¬ 
sive relays are used to connect two bodies in interaction. 
In this paper, we study the thermodynamic performance 
of such a 3B system and demonstrate the strong potential 
of these near-held heat engines for energy harvesting. 


The properties of the thermal radiation driving an 
energy-conversion process depend on the distribution and 
the number of bodies interacting with the converter. Let 
us compare the operating modes of two-body (2B) and 
3B radiative heat engines, which are both sketched in 
Fig. [TJ In a 2B heat engine a hot source at temperature 
Th radiates towards a converter which is assumed to be in 



FIG. 1. Sketch of a heat engine with a hot source at tem¬ 
perature Th and a cold sink at temperature T c < Th that 
provides a usable work flux W by converting near-field ther¬ 
mal radiation energy. The cold sink receives a heat flux Q. 
(a) Three-body: One of the bodies (emitters) is thermalized 
with the source and another one with the sink, while a passive 
intermediate body (of width 5) is placed between them. The 
net energy and entropy fluxes on the cold body are < 1 > 3 b and 
4/3B, respectively, (b) Two-body: The intermediate body is 
removed. The net energy and entropy fluxes on the cold body 
are <I> 2 b and \1/2B, respectively. The distance d between the 
bodies is indicated in both cases. 


contact with a cold sink at temperature T c < Th- In the 
3B configuration, a passive intermediate body of width 
8 is placed between the source and the sink. Note that 
this passive relay is maintained at the same separation 
distance d from both the source and the sink as the cav¬ 
ity width in the 2B system. Hence, we do not introduce 
in the 3B heat engine an exaltation mechanism which re¬ 
sults from a simple reduction of distances. Moreover, we 
assume that the intermediate body reaches a local equi¬ 
librium temperature T r , an assumption that is justified 
in practical applications for the sizes of the body that 
we consider here. This temperature T r is not arbitrary. 
It is taken such that the net energy flux that the inter¬ 
mediate body exchanges with the hot and cold bodies 
vanishes. Hence, T r is an implicit function of Th, T c , and 
the two parameters d and 8 that specify the geometry 
of the problem. As a consequence, the energy flux radi¬ 
ated by the hot body coincides with the flux received by 
the cold body. This ensures that the energy supplied to 
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the system comes only from the hot source, since, under 
these conditions, a thermostat at T r in contact with the 
passive relay will provide a vanishing net energy flux. 

The planar 2B and 3B structures considered here have 
an infinite transversal extension. When all the bodies 
are separated by vacuum, the net energy flux on the cold 
body can be written as an integral over monochromatic 
contributions of frequency w, which in the near-field 
regime is given by = f 0 °° d, 5) (i = 2,3) 

with 


</>3b(w, d, S) = hw^2 

•J CK 


d 2 k 
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for the 3B configuration |15| and 
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in the 2B case 0|- Above, we have intro¬ 
duced n a p(oj) = n a (ui) — np(uj), where n a {u>) = 
^ e htj/k B T a _ are distributions of photons at 
equilibrium temperature T a with a = h,r, c, being 
Boltzmann’s constant and 2irh Planck’s constant. In 
Eqs. © and m, the sum runs over the polarizations 
j = s,p, and the integration is carried out over the com¬ 
ponents of the transverse wave vector k. = ( k x ,k y ) with 
k = |k| > aj/c, c being the speed of light in vacuum. 
This means that only the dominant near-field contribu¬ 
tion of evanescent waves is taken into account, whereas 
the contribution of the propagating modes with k < w/c 
is neglected. The transmission coefficients of the 2B and 
3B systems are defined as mm 

Tj hr) = 4 | tJ | 2 Im (pf Im (pf e~ A ^ d / \Df''Df | 2 , 

T} rc) = 41m (pf) Im (pf) e ~ 2 ^ d / \Df c \ 2 , (3) 

7} hc) = 41m (pf Im (pf e ~ 2 ^ d /\Df\ 2 , 

where k z = fu 2 /c 2 — k 2 is the normal component of 
the wave vector while Df = 1 — p b pje 2lkzd , Df c = 
1 - pfp'-e 2lk ‘ d and Df = 1 - p l ]p c j e 2ik * d are the Fabry- 
Perot-like denominators. Here p b j = pfu,K), p £ = 
Pj(ui,K,d), and pj = p c fu>,n) are the reflection coeffi¬ 
cients of the hot, intermediate, and cold bodies, respec¬ 
tively, tJ = tJ(ui,k, 5) are the transmission coefficients of 

the intermediate body, and pf = P^j+ff 2 p b e 2q,kzd / Df 
are the reflection coefficients of the hot and the interme¬ 
diate bodies considered as a single entity. 

The expressions © for the transmission coefficients 
show that the three bodies are coupled together due to 
multiple interaction mechanisms resulting in their non¬ 
trivial optical properties. Moreover, in view of Eqs. © 
and ©, the energy flux <fq B can be written as a sum 


FIG. 2. (a) Ratio of the maximum work flux in the 3B con¬ 
figuration to the maximum work flux in the 2B configuration, 
VV 3 B/W 2 B, as a function of the separation d and width of the 
intermediate body S for Th = 400 K and T c = 300 K. A region 
of amplification due to 3B photon tunneling is clearly appre¬ 
ciated. (b) Efficiency ratio J 73 b/?? 2 b in the same conditions. 
In (c) and (d), the plots show the corresponding maximum 
work flux and efficiency in the 3B configuration. 


over the contributions stemming from the different con¬ 
stituents of the system which are in local thermal equi¬ 
librium, i.e. < ^’a B ' 1 (T a ) with a = h,c for 2B 

and a = h,r,c in the 3B case. The fluxes 4>a B ^(T a ) de¬ 
pend only on the local equilibrium temperatures T a of the 
constituents through the distribution functions n a (u >). 
Hence, the partial entropy fluxes 'Pa iB ' l (T , a ) carried by 
the thermal fields generated by the different constituents 
are given by [l3| 


*£ B) (T a ) = j V1 A$£ B ) (r '). (4 ) 

Therefore, the net entropy flux on the cold body reads 

4' 1 b = E^ B) ( T “) = l 2^iB(w,d,*), (5) 


where the spectral entropy fluxes take the form 

ip 3B (u,d,5) = k B J2 I (6) 

j J CK >UJ ) 

x m hr (u})Tj hr \uj,K,d,6) + rn rc (u})Tf c \u},K,d,d) , 
V>2B(w,d) = fc B ^ [ -£^m hc (u)Tj hc) (u,K,d), (7) 
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with m a p(oj ) = m a (ui) — rrip(uj) and m Q (w) = 

[1 + n a (ui)] In [1 + n a (w)] — n a (cu) lnn a (w). These last 
two relations are strictly valid only if the temperature 
dependence of the material properties of the constituents 
of the 2B or 3B system can be neglected in the considered 
range of working temperatures. 

Once energy and entropy fluxes are known, the ther¬ 
modynamics of the energy-conversion process can be an¬ 
alyzed as follows (the 3B and the 2B configurations will 
be discussed simultaneously). First of all, notice that, 
due to the difference of temperatures between the bod¬ 
ies, the transport of heat through the cavity proceeds 
irreversibly and entropy is generated at a certain rate, 
say, 'F g . This entropy production \k ff accounts for dissi¬ 
pative processes in the thermalization of excited electrons 
at the surface of the cold body 22]. Since the bodies are 
thermalized, in particular, a heat flux Q is transferred 
isothermally to the cold sink; we assume that this trans¬ 
ference is done reversibly and, thus, Q = T c (dbe + 'kg). 
In this scheme the heat engine can be considered as en- 
doreversible, as discussed in Q for the 2B problem in the 
far field. Taking into account the balances of energy and 
entropy fluxes, the work flux that can be delivered by the 
engine reads W = — T c (dbe + 'kg). Since *kg > 0, 

the maximum work flux or thermodynamic availability is 
given by Wm = ^iB — T c ^i B . In addition, considering 
<kiB as the input energy flux, the efficiency of the engine 
is given by r]i B = W’/^iB- According to this, an upper 
bound for the efficiency can be obtained by computing 
the ratio fj lB = 

In the configuration we analyze, the hot and cold 
bodies are two 5-/xm-thick silicon carbide (SiC) sam¬ 
ples [23] that support a SPP with a resonance at wspp — 
1.79 x 10 14 rad/s. Moreover, as in Ref. jUj, we use in 
the 3B configuration for the intermediate slab a metal¬ 
like medium which supports a surface mode (a plasmon) 
at the same frequency wspp. Figure [2](a) shows the ratio 
of the maximum work flux in the 3B system to the max¬ 
imum work flux in the 2B system, i.e., with and without 
an intermediate relay. It can be seen that a 3B engine 
can produce about 60% more work than a classical 2B 
system. If the width of the intermediate body becomes 
sufficiently large, the 3B interaction disappears (in the 
near-field regime), and both cavities, located between the 
source and the intermediate relay and between the relay 
and the sink, become independent. Then, the work pro¬ 
duction by the 3B heat engine becomes comparable to 
or even smaller than the one of a 2B engine. As for the 
efficiency of those engines, we see in Fig. [2](b) that they 
are comparable in both configurations provided the sepa¬ 
ration distances are large enough compared to the width 
of the intermediate slab. It is interesting to note that the 
2B efficiency seems always to be larger than the 3B effi¬ 
ciency even in the parameter range where the extracted 
work of the 3B system exceeds that of the 2B system. 

In order to get some insight on these results, we plot 
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FIG. 3. (a) Transmission coefficients weighted by the pho¬ 
ton distributions taking d = 500 nm and S = 667 nm. We 
plot = 10 22 x ( nhrTp hr) + n rc Tp rc) j for Th = 400K 

and T c = 300K, for which T r = 357.01 K. In (b), /( lo,k) is 
shown for Th = 200 K and T c = 100 K with T r = 180.54 K. 
(c) Spectral energy fluxes 03 b and (p 2 B in 3B and 2B configu¬ 
rations, respectively, corresponding to the same setting used 
in (a), (d) Spectral energy fluxes 4>3b and 4> 2 B corresponding 
to (b); the inset shows the same spectra in log scale. In (e) 
and (f) we plot the spectral entropy flux (multiplied by the 
temperature of the sink) T c ips b and the spectral work flux 
W 3 B = </> 3 B — T c i/j 3 b corresponding to the same setting used 
in (a)-(c) and in (b)-(d), respectively. 


in Figs. [3]]a) and Elb) the transmission coefficients for 
p-polarized waves (the main contribution) in the (w, k ) 
plane associated to the 3B engine by weighting them with 
the corresponding photon distribution functions. As a 
first observation, we see the presence of different surface- 
mode branches of the four coupled surface modes (sym¬ 
metric and antisymmetric modes) in the 3B system [24| 
around the surface-mode resonance frequency of SiC. For 
these surface modes, the transmission is apparently high. 
These branches support high transmissions for large wave 
vectors, which means that a large number of modes con¬ 
tribute to the heat transfer in this spectral region W- 
The closer the frequency of the surface mode gets to 
Wien’s frequency of the heat source uw = 2.82kBTh/h, 
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the higher the number of excitations that contribute to 
the transfer. Accordingly, if the hot body is cooled down 
to a temperature for which is far from wgpp, the 
modes in the region around the SPP stop to contribute 
to the transfer, as can be seen in Fig. [3|b). The spec¬ 
tral energy fluxes plotted in Fig. ©c) and ©d) corrob¬ 
orate this tendency. Since the 3B photon-tunneling en¬ 
hancement occurs in the SPP region, we thus observe in 
Fig. Sc) an increase of the quasimonochromatic spec¬ 
tral energy flux <f> 3 B as compared with Further¬ 

more, the spectral entropy flux and the spectral 
work flux W 3 B = 4> 3B ^ T c ip 3B are also peaked around 
the SPP frequency. As shown in Fig. Se), the neg¬ 
ative entropic term drastically reduces the monochro¬ 
matic contribution w^b to the thermodynamic availabil¬ 
ity W 3 B(d, <5) = / 0 °° ^W 3 b(w, d, S). This entropic term 
represents a non-negligible energy flux that the system 
transfers to the cold sink, thus diminishing the amount 
of usable work production. 

Finally, we study the maximal work flux W,b that can 
be extracted from such a 3B heat engine compared with 
that of a 2B heat engine. To this end, we fix T c = 300 K, 
while the temperature of the heat source Th is varied. 
For the thickness of the vacuum gaps and the intermedi¬ 
ate passive relay we choose d = 100 nm and 5 = 133 nm 
for which W 3 B/W 2 B is maximum when Th has reached a 
temperature of 400 K, as shown in Fig. ©a) (in that case, 
the ratio W 3 B/FV 2 B slowly increases for increasing Th). 
The results are plotted in Fig. ©a). They show that the 
discrepancy between W 3 B and W 2 B grows monotonically 
with respect to temperature while the 2B and 3B effi¬ 
ciencies remain very close to each other (t^b > V 3B ). In 
addition, the dependence of the work fluxes and efficien¬ 
cies on the distance d for the optimal S are presented in 
Fig. ©b), where the temperatures are set to Th = 400 K 
and T c = 300 K. This example illustrates that in a 3B 
system the energy flux and the maximal work flux that 
can be extracted are enhanced by the interactions of the 
surface modes in the hot and cold body with that of the 
intermediate relay. 


The maximum transfer in a 3B configuration takes 
place when the transmission coefficients attain their max¬ 
imum value. The theoretical limit is thus achieved by 

as also occurs for 


the condition T- hr ' > = 1= 1 


2B systems [18|, |20|, |2l| when T^ hc ^ = 1. Using this 


in ©> m, ©^ and m and taking into account a cut¬ 
off wave vector k c ,iB w/c, for which the modes are 

effectively confined, we get 4>^j ax = (T 2 — T 2 ) and 

= 2&b (: T h -T c ), where & B = < iB fc|/24A No¬ 
tice that maximizing the energy flux implies also that 
the flow of entropy per channel is maximum [20| . The 
maximum work flux is thus kV I 1 § ax = £ib(T)i — T c ) 2 , and 
in consequence the upper bound for the efficiency reads 
V 2 B X = Cb x = (T h - T c )/(T h + T c ). Therefore, the effi¬ 
ciencies for the 2B and 3B systems are equal. However, 



FIG. 4. (a) Maximum work fluxes W 3 B and W 2 B for the 

3B and 2B configurations, respectively, as a function of the 
temperature of the hot source Th- The temperature of the 
cold sink is set to T c = 300 K, the separation distance to 
d = 100 nm, and the width of the intermediate body to <5 = 
133 nm. The inset shows the corresponding upper bounds for 
the efficiency, fj 33 and f) 2 b ■ (b) Dependence of the maximum 
work fluxes and efficiencies on the separation distance d for 
fixed temperatures and the optimal S. 


we remark that the difference between a 3B and a 2B sys¬ 
tem is manifested through K Cj iB- The cutoff wave vector 
in a 3B system can be larger than that of the 2B config¬ 
uration as shown in our numerical examples. Although 
the efficiencies and the ratios are the same 

in the 3B and 2B system, it follows that 


W max (Rmax 

3B _ _j_3B 

W max ^max 

2B ^2B 


K c,3B 

K c,2B 


> 1 . 


( 8 ) 


Hence, a larger maximum work flux in the 3B system is 
due to the larger energy flux which, in turn, results from 
the larger number of contributing modes. 

For the general case of an A-body (AB) engine, we 
can always write the heat flux exchanged between the 
different parts of the system by using the Landauer for¬ 
malism derived in Refs. 18320!], the transmission co¬ 


efficients being simply related to the scattering of the 
electromagnetic field radiated by each part of the sys¬ 
tem. In this case, a suitable cutoff k C: nb has to be 
considered. When k c ,nb > k c , 2 B for A > 2, an en¬ 
hancement in the maximum work flux W™g x is thus ex¬ 
pected with respect to the 2B case, while the efficiency 
Vnb? = (Th ^ T c )/(Th + T c ) remains the same for all N. 
Furthermore, the heat conductance in this N B endore- 
versible engine is, in general, nonlinear in the difference of 
temperatures AT = Th — T c . However, the linear regime 
can be achieved if AT <C Th, and, under this assump¬ 
tion, the efficiency becomes fj ™|f = AT/(2Th), which is 
half the Carnot efficiency. For engines with linear heat 
conductance such as the Novicov engine, the efficiency is 
given by t/ca = 1 - \fT c /T h (H), which for AT < T h 
becomes tjca = AT/(2T^). We thus observe that the ef¬ 
ficiency of our engine, in the appropriate limit, coincides 
with that of a Novicov engine. 

In conclusion, we have demonstrated that the thermo- 
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dynamic performance of 3B near-field heat engines can 
substantially overcome that of 2B systems. Our results 
pave the way for a generation of nanoscale energy con¬ 
verters driven by the physics of many-body interactions 
instead of the conventional two-body interactions. In ad¬ 
dition, this work provides perspectives for investigating 
the thermodynamics of systems with long-range electro¬ 
magnetic interactions. 
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